A GRADIENT FLOW APPROACH TO A THIN FILM APPROXIMATION 

OF THE MUSKAT PROBLEM 



PHILIPPE LAURENgOT AND BOGDAN-VASILE MATIOC 



Abstract. A fully coupled system of two second-order parabolic degenerate equations arising as a thin 
film approximation to the Muskat problem is interpreted as a gradient flow for the 2-Wasserstein distance 
in the space of probability measures with flnite second moment. A variational scheme is then set up and 
is the starting point of the construction of weak solutions. The availability of two Liapunov functionals 
turns out to be a central tool to obtain the needed regularity to identify the Euler-Lagrange equation in 
the variational scheme. 



1. Introduction 

The Muskat model is a free boundary problem describing the motion of two immiscible fluids with 
different densities and viscosities in a porous medium (such as intrusion of water into oil). Assuming 
that the thickness of the two fluid layers is small, a thin film approximation to the Muskat problem 
has been recently derived in [10] for the space and time evolution of the thickness / = f{t,x) > and 
g = g{t, x) > of the two fluids {f + g being then the total height of the layer) and reads 

f dtf = il + R)dAfdJ) + RdAfd.g), ,,,,,, ™ 

{ {t,x) e (0,oo) X M, (1.1a) 

[ dtg = Rf^dx (gdxf) + i?^«9x (gdxg) , 

supplemented with the initial conditions 

/(0) = /o, 9(0) = go, x£R. (1.1b) 

Here, R and R^ are two positive real numbers depending on the densities and the viscosities of the fluids. 
Since / and g may vanish, (1.1a) is a strongly coupled degenerate parabolic system with a full diffusion 
matrix due to the terms dx{fdxg) and dx{gdxf)- There is however an underlying structure which results 
in the availability of an energy functional 

£U,9):=\ I [f+Rif + gf] dx, (1.2) 

which decreases along the flow. More precisely, a formal computation reveals that 

-£{f,g) = - f {{1 + R)dxf + Rdxgf + RR^, g {dxf + dxgf dx . (1.3) 
at j]^ I J 

A similar property is actually valid when (1.1a) is set on a bounded interval (0, L) with homogeneous 
Neumann boundary conditions: in that setting, the stationary solutions are constants and the principle 
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of linearized stability is used in [10] to construct global classical solutions which stay in a small neighbour- 
hood of positive constant stationary states. Local existence and uniqueness of classical solutions (with 
positive components) are also established in [10] by using the general theory for nonlinear parabolic 
systems developed in [2]. Weak solutions have been subsequently constructed in [9] by a compactness 
method: the first step is to study a regularized system in which the cross-diffusion terms are "weakened" 
and to show that it has global strong solutions, the proof combining the theory from [2] for the local 
well-posedness and suitable estimates for the global existence. Some of these estimates turn out to be 
independent of the regularisation parameter and provide sufficient information to pass to the limit as 
the regularisation parameter goes to zero and obtain a weak solution to (1.1a) in a second step. A key 
argument in the analysis of [9] was to notice that there is another Liapunov functional for (1.1a) given 

by 



which evolves along the flow as follows: 



R 

/In / + — ghig 



dx . 



(1.4) 



d_ 



n{f,g) = -f [m\^ + R\d^f + d.,g\ 



dx . 



The basic idea behind the above computation is to notice that an alternative formulation of (1.1a) is 




d,.[f d^{{l + R)f + Rg)], 
Rfi dx [g dx if + g)] , 



{t,x) G (0,oo) X 



so that it is rather natural to multiply the /-equation by In / and the (^-equation by Ing and find nice 
cancellations after integrating by parts. In this note, we go one step further and observe that a concise 
formulation of (1.1a) is actually 



dtf = dx 



R_ 

Rii 



dtg = dx 



f (^0(/'5) 
g dx 



{t,x) G (0,oo) X M, 



:i.5) 



which is strongly reminiscent of the interpretation of second-order parabolic equations as gradient flows 
with respect to the 2-Wasserstein distance, see [3, Chapter 11] and [18, Chapter 8]. Indeed, since 
the pioneering works [12] on the linear Fokker-Planck equation and [15, 16] on the porous medium 
equation, several equations have been interpreted as gradient flows with respect to some Wasserstein 
metrics, including doubly degenerate parabolic equations Jl], a model for type-II semiconductors [4], 
the Smoluchowski-Poisson equation [5], some kinetic equations [6,8], and some fourth-order degenerate 
parabolic equations [14], to give a few examples, see also [3] for a general approach. As far as we know, 
the system (1.5) seems to be the first example of a system of parabolic partial differential equations which 
can be interpreted as a gradient flow for Wasserstein metrics. Let us however mention that the parabolic- 
parabolic Keller-Segel system arising in the modeling of chemotaxis has a mixed Wasserstein-L^ gradient 
flow structure [7]. 

The purpose of this note is then to show that the heuristic argument outlined previously can be made 
rigorous and to construct weak solutions to (1-1) by this approach. More precisely, let /C be the convex 



subset of the Banach space L^(M, (1 + x'^)dx) n L^(M) defined by 

/C := |/i G (1 + x'^)dx) n L2(E) : /i > a.e. and y = l| , (1.6) 

and consider initial data {fo,9o) G /C2 := /C x /C. We next denote the set of Borel probabihty measures 
on M with finite second moment by 7-'2(IR) and the 2— Wasserstein distance on V2(^) by Recall that, 
given two Borel probability measures fi and z/ in 7^2(1^)! 



Wi{fi,i^):= inf / \x -y\'^d7r{x,y) 



where n(/i, 1^) is the set of all probability measures vr G V{M?) which have marginals fi and u, that is 
■7r[A X M] = n[A] and 7r[M x B] = for all measurable subsets A and i3 of M. Alternatively, vr G n(;U, i/) 
is equivalent to 

{cl){x) + ?/^(y)) ?/)=/" (/.(x) + /" duiy) for all V) G l1(M; R^). 

With these notation, our result reads: 

Theorem 1.1. Assume that R > 0, > 0. Given r > and {fo,go) G IC2, the sequence (/",(7")„>o 
obtained recursively by setting 

{flg^,) := (/o,<7o), (1.7) 



with 



5."+^) := inf -^(n,t;), (1.8) 

^(n,^;):=^ (^l{u,n + ^^Wi{v,a^,)^+£{u,v), (u,^;)g/C2, (1.9) 

is well-defined. Introducing the interpolation [f^^g^) defined by 

fr{t) := and gr{t) := g"^ for t G [nr, {n + 1)t) and n > 0, (1.10) 

there exist a sequence {Tk)k>i of positive real numbers, Tk \ 0, and functions {f,g) : [0, 00) — IC2 such 
that 

ifr^gr,) {f,g) inL^{{0,T)xR;R^)forallT>0. (1.11) 

Moreover, 

{i) {f,g) G L-(0,oo;L2(M;R2)), (a,/,^,^) G L\0,t; H^{R-R^)), 
(ii) {f,g) G C([0,oo);ii-3(M;R2)) i^zi/, (/,5)(0) = (/o,5o), 
and the pair {f,g) is a weak solution of (1.1) in the sense that 

ft 



I fit) ^dx- [ foCdx+ [ [ f{a) [(1 + R)d,f + Rd^g] {a)d,^ dxda = 0, 

JR JR JO JR 

[ g{t)idx- f go^dx + Rf, [ [ g{a) [d^f + d^g) {a)dx£,dx da = Q , 

V JR JR Jo JR 



(1.12) 
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for all 1^ G C^(M) and t > 0. In addition, {f,g) satisfy the following estimates 

(a) n{f{T),g{T))+ r f [\dj\^ + R\d^{f + g)\''] dx dt <U{fo,go) , 

Jo Jm. 

(b) £{f{T),g{T)) + \ f I \f ((1 + R)d.,f + Rd^^gf + RR^^g^d^f + d^gf] dxdt < £ifo,go) , 

^ Jo Jm. ^ 

for a.e. T £ (0,oo), £ andU being the functional defined by (1.2) and (1.4), respectively. 

Let us briefly outline the proof of Theorem 1.1: in the next section, we study the variational problem 
(1.8) and the properties of its minimizers. A key argument here is to note that the availability of 
the Liapunov functional (1.4) allows us to apply an argument from [14] which guarantees that the 
minimizers are not only in L^(]R;]R^) but also in i7-^(M;M^). This property is crucial in order to derive 
the Euler-Lagrange equation in Section 2.2. The latter is then used to obtain additional regularity on the 
minimizers, adapting an argument from [15]. Convergence of the variational approximation is established 
in Section 3. Finally, three technical results are collected in the Appendix. 

As a final comment, let us point out that we have assumed for simplicity that the initial data /q and go 
are probability measures but that the case of initial data having different masses may be handled in the 
same way after a suitable rescaling: more precisely, let (/o,ffo) S L^(M) H L^(R, (1 + x'^)dx) and denote 
a solution to (1.1) by (/,g). Setting F := //||/o||i and G := g/\\go\\i and recalling that ||/(t)||i = ||/o||i 
and ||5'(i)||i = Hffolli for all i > 0, we realize that {F,G) solves 



^ dtF = d,[F d,{{l + R)7]^F + RG)] 



\\9o\\i 
R 



{t, x) G (0, oo) X 



^mII/oI 



1 



dtG = d, [G d, {RF + Rrj-^G)] , 



with := ||/o||i/|bo||i and initial data (Fo,Go) := (/o/||/o||i,5o/||5o||i) G ^2- The corresponding 
variational scheme then involves the functional 

1^ ( tAt Wi{u,Fo) + -^W^{v,Go))+^ Ml + ^ + (n,^)G/C2, 
to which the analysis performed below (with r] = 1) also applies. 

2. A VARIATIONAL SCHEME 
Given r > and (/o,9o) G ^2, we introduce the functional 



jr^{u,v):=—i^W^{u,fo) + —W^{v,go)]+£{u,v), {u,v)e)C2, (2.1) 
and consider the minimization problem 



inf Tr{u,v). (2.2) 

{u,v)GK.2 



2.1. Existence and properties of minimizers. Let us start by proving that, for each (/o,(7o) £ 
the minimization problem (2.2) has a unique solution in /C2- 

Lemma 2.1. Given (/o,(7o) £ ^2 a^^^ r > 0, there exists a unique minimizer {f,g) G /C2 0/ (2.2). 
Additionally, {f,g) G F1(M;]R2) ^y^^/j 



\dJ\\l + R\\d^{f + g)\\l<- 

T 



(2.3) 



where 



H{h) := / hln{h) dx for h G L^(]R) smc/i i/iai /i > a.e. and hln{h) G L^(M). (2.4) 
Jr 

Recall that, if /i G /C, then hlnh G L^(M) (see Lemma A.l below) so that the right-hand side of (2.3) 
is well-defined. 

Proof. The uniqueness of the minimizer follows from the convexity of IC2 and W2 and the strict convexity 
of the energy functional £. 

We next prove the existence of a minimizer. To this end, pick a minimizing sequence (uk, Vk)k>i G IC2- 
There exists a constant C > such that 

Ikfclb + Ibfclb < C, k>l, (2.5) 
W2{uk,fo) + W2{vk,go) < C, k>l. (2.6) 

From (2.5) we obtain at once that there exist (/, g) G L^(M; M?) and a subsequence of {uk,Vk)k>i (denoted 
again by (iifc,ffc)fc>i) such that 

Uk^f and Vk^g inL2(M). (2.7) 

Let us first check that (/, g) G 1C2- Lideed, the nonnegativity of / and g readily follows from that of 
Uk and Vk by (2.7) while integrating the inequality x"^ < 2y^ -|- 2\x — with respect to an arbitrary 
vr G n(iifc,/o) yields 

/ Uk{x)x'^ dx = j x'^d7r{x,y) < 2 y'^ d7r{x,y) + 2 / \x — y\'^ dTr{x,y) 
Jr Jr^ Ji82 J^2 

< 2/ fo{y)y^dy + 2 [ \x - yl"" dTT{x,y) , 

which implies by virtue of (2.6) that 

/ Uk{x)x'^dx<2 / fo{x)x'^dx + 2W^{uk,fo) <C , k>l. (2.8) 
Jr Jr 

Similarly, 

Vk{x)x^dx<C, k>l. (2.9) 



Owing to (2.5), (2.8), and (2.9), we deduce from the Dunford-Pettis theorem that {uk)k>i and {vk)k>i 
are weakly sequentially compact in L^(R). We may thus assume (after possibly extracting a further 
subsequence) that ^ f and f ^ ^ (7 in L^(]R), whence 

/ f{x)dx= lim / U}^{x)dx = l and / g{x) dx = lim / Vk{x)dx = l. 
Jr k^'^jR Jr k^°°jR 
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Finally, combining (2.5), (2.8), and (2.9) with a truncation argument ensure that / and g both belong 
to L^(]R, (1 + x'^)dx). Summarising, we have shown that {f,g) G JC2- 
The next step is to prove that 

^r{f,g) = , inf Fr{u,v) . 
Indeed, on the one hand, the weak convergence (2.7) implies that 

£{f,9) < liuiiniS {uk,Vk). 

On the other hand, we recall that the 2-Wasserstein metric W2 is lower semicontinuous with respect to 
the narrow convergence of probability measures in each of its arguments, see [3, Proposition 7.1.3], and 
the weak convergence of {uk,Vk)k>i in L^(]R;M^) ensures that 



Consequently, 



Wiif, fo) < lim inf Wiiuk, fo) and W^ig, go) < lim inf WUvk,go) ■ 

fe— ^-oo fc— >oo 



J^T{f,g) < liminf JV(nfc,-Ufc) with {f,g) G IC2 



SO that {f,g) is a minimizer of Jv in IC2- 

As a final step, we show that / and g belong to H^{W). To this end, we follow the approach developed 
in [14] and take advantage of the availability of another Liapunov function as already discussed in the 
Introduction. More precisely, denote the heat semigroup by {Gt)t>Oy that is, 

{Gth){x) := ^= [ expf- '^"^'' ) h{y) dy , (t, x) G [0, 00) x M , 

for h G L^(IR). Since (/, 5) G JC2, classical properties of the heat semigroup ensure that {Gtf,Gtg) G JC2 
for all t > 0. Consequently, J-'r{f,g) < J-r{Gtf, Gtg) and we deduce that 



£{f,g)-g{Gtf,Gtg)<^ 



[WiiGtfJo) - WiifJo)) + ^ {WiiGtg,go) - Wi{g,go)) 



(2.10) 



£{Gtf, Gtg) = [ [Gtf dtGtf + R [Gtf + Gtg) dt{Gtf + Gtg)] dx = -\\d,Gtf\\l - R \\d^Gt{f + g)\\l 
Jr 



for all t > 0. Moreover, for all t > 0, we have 
d 
dt 

and by integration with respect to time we find that 

7 f[\\d.Gsfg + R\\d.Gs{f + g)f2] rf,< ^(/^g)-^(gt/,Gtg) foralH>0. 
t Jo t 

Since s 1— )■ ||(?2,.Gs/i||2 is non-increasing for h G L^(M) we end up with 

Wd.GtfWl + R\\dMf + g)\\l< ^^f^9)-nGtf,Gtg) ^ ^ ^ ^2.11) 

We recall now some properties of the heat flow in connection with the 2-Wasserstein distance W2, 
see [3,8, 16], these properties being actually collected in [14, Theorem 2.4]. The heat flow is the gradient 
flow of the entropy functional H given by (2.4) for 14^2 and, for all (h, h) G IC2, we have [3, Theorem 11.1.4] 

-—W^{Gth, h) + H{Gth) < H{h) for a.e. t > 0. (2.12) 



Choosing {h,h) = (/, /o) and {h,h) = {g,go) in (2.12), we obtain 



2di 



Wi{Gtf,fo) + ^Wi{Gtg,go) 



R 



< Hifo) - H{Gtf) + — {H{go) - H{Gtg)) 

rill 



for a.e. t > 0. Integrating the above inequahty with respect to time and using the time monotonicity of 
s I— 7- H{Gsf) and s i— ?• H{Gsg) give 



1 
2 

< 
< t 



WiiGtfJo) - W^ifJo) + 4- (^l(Gt5,5o) - W^ig,go)) 



Hifo) - H{Gsf) + 



R_ 

Rii 



[Higo) - H{Gsg)) 



ds 



R 



H{fo) - H{Gtf) + — {H{go) - H{Gtg)) 



Gathering (2.10), (2.11), and (2.13), we find 



\\dccGtf\\l + R\\d,Gt{f + , 



< 



Hifo) - H{Gtf) + 4- (Higo) - H{Gtg)) 
rill 



(2.13) 



(2.14) 



for i > 0. As a direct consequence of (2.14) and the boundedness from below (A. 2) of H in /C, {dxGtf)t>o 
and {dxGtg)t>o are bounded in L^(]R) and converge to dxf and dxg^ respectively, in the sense of distri- 
butions as t — )• 0. This implies that both / and g belongs to H^[W) and we can pass to the limit as t — t- 
in (2.14) to obtain the desired estimate (2.3) and finish the proof. □ 

2.2. The Euler-Lagrange equation. We now identify the Euler-Lagrange equation corresponding to 
the minimization problem (2.2). 

Lemma 2.2. Given {fo,go) G ^2 CLud r > 0, the minimizer {f,g) of Tr in IC2 satisfies 



- U{f-k)dx+ f [((1 + R) f dxf + Rf dxg) dxi] dx 
- I C {g - go)dx + Rf, I [{gdxf + gdxg) dxi] dx 



< 



< 



2t 
2r 



-wlifjo) 

- Wl{g,go) 



(2.15) 
(2.16) 



forieC^m. 

Proof. To derive (2.15)-(2.16) we follow the general strategy outlined in [18, Chapter 8]. According to 
Brenier's theorem [18, Theorem 2.12], there are two convex functions : M — )■ M and ^ : R — R which 
are uniquely determined up to an additive constant such that 



\x - dxip{x)\'^ fo{x) dx 



inf 

T#fo=f 



where the infimum is taken over all measurable functions T 
i.e. satisfying 



\x-T{x)\^ fo{x)dx, (2.17a) 
pushing /o forward to / (/ = T#/o), 



f{x) dx 



B 



and 



Wi{g,go) 



IT-^{B) 

\x - dxip{x)\'^ go{x) dx 



foix) dx for all Borel sets B of R, 

X — S{x)\'^ go{x) dx. 



inf 

S#go=9 



(2.17b) 
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We pick now two test functions ij and ^ in C^(M) and define 

/, := ((id +eO o d^ip)#fo = (id +eO#f and g, := ((id +erj) o a,V)#^o = (id +er])#g (2.18) 

for each e G [0, 1], where id is the identity function on M. To ease notation we set 

re:=id+e^ and Se:=id+er], (2.19) 

and observe that there is eo smaU enough (depending on both ^ and rj) such that, for e £ [0,eo], and 
Sg are C°°— diffeomorphisms in M. Then, by (2.18), we find the identities 

fe= JrX-^ '^ = 7ri^' ^e(0,eo]. (2.20) 

dxTe o Te OxSe o bg 

Observing that ||/e||i = ||/||i = \\gs\\i = \\g\\i = 1 and 

we clearly have {fe,ge) G ^2 for all e G (0,eo] a-nd thus J-rifjg) ^ -^r(/e)ffe)- Consequently, 
1 



< 
- 2r 



Wiife, /o) - Wiif, fo) + -|- {Wiige, go) - W^{g, go)) 



+ £{fs,gs)-£{f,g). (2.22) 



Concerning the energy £, it follows from (2.21) that 

2{£ifs,gs) - £{f,g)) = {1 + R) I! + RII + 2R If , (2.23) 

with 

Jr \dxTe{x) ) Jr \d:^Se{x) J 



/|:= / ifege-f g){x)dx. 

We now consider the three integrals in the right-hand side of the relation (2.23) separately: since 



/f = -£ / ±^ dx and n = -e / g^ dx , 

it readily follows from Lebesgue's dominated convergence theorem that 

lim ^ = - / d^i f dx and lim ^ = - / 9^r/ g'^ dx . (2.24) 
We next turn to the term /g involving both / and g and split it in two terms 2/3 = /g^ + with 
III ■■= I ife + f) i9e - 9) dx and /fa := / {g, + g) (/, - f) dx . 

By (2.20), 

Ik=f(^^^f) (^^-9) dx 

7r V d^T, oTe^ V V o Se^ V 

= [ ( J°X^l\ +foSs) {g - {g o Se) d,S,) dx 

= [ +foSe){g-goSe)dx-e [ d.r^ {g o S,) f -1^10^ + f o S,) dx. 

jR\dxTeOTe o be J Jr \dxTe o Te o be J 



On the one hand, invokmg Lemma A. 2 (with {h,() = {g,il))j '^^ know that {g — g o Sg)/e —rfd^g 
in L^(]R) as e — 7- 0. On the other hand, using again Lemma A. 2 as weU as Lemma A. 3, we have that 
f o Ss — > f and / o o — > f in L^(]R) as e — )• 0, and so does / o o Ss/{dxTir o Tf^ o Se) owing 
to the uniform convergence of (9x^e)e to 1. Consequently, 



lim:|i=-2 / fdx{rjg)dx, 



and similarly 

lim-^=-2 /" gd.mdx. 
Gathering (2.24)-(2.26), it follows from (2.23) that 



(2.25) 



(2.26) 



lim 

e-i>0 



£Ue,ge)-£if,g) 



{l + R) I ^d^^dx-R I '-d^iidx-R I [f dx{iig) + g dx{U)] dx 



(2.28) 



(2.29) 



(2.27) 

To handle the terms of (2.22) involving the Wasserstein distance, we argue as in [18, Section 8.4] and 
write 

W^{feJo)< I \id-TeodM'' kdx= I \ id-dx^-eiodM'' hdx 
Jr Jr 

= / \id-dx^p\'^ fodx -2e l iid-d^^) {^o d^^) f^dx + / \^odx^\^ f^dx, 
Jr Jr Jr 

from which we deduce, according to the definition of dx^p, 

Wiif,, /o) < Wiif, /o) - 2e / (id -dx^) (^ o dx^) fo dx + e^ [ |^ o dx^f /o dx, 

Jr Jr 

and similarly 

Wi{ge,go) < W2{g, go) - 2e / (id -dxi^) (r? o dx^p) godx + e^ / jr? o dx^'^ go dx. 

Jr Jr 

Summing (2.27), (2.28), and (2.29), we obtain by dividing (2.22) by e and letting e ^0 that 

- [ (id -dx^p) (^ o dxif) fodx + ^ [ (id -dxi') (jj o dxi^) go dx 
T Ur -"-m -'k 

+ {1 + R) [ dx^^dx + R [ dxr]^dx + R [ [f dx{r,g) + g dxiU)] dx < 0. 
^ Jr ^ Jr 

Since the relation is valid for (^,7/) as well as for (— ^, —r]), we end up with 

- [ (id-dxif) {iodx^p) fodx + ^ [ (id -5^^) {rj o dxi)) godx 
T Ur -"-m Jr 

+ {l + R) [ dxC^dx + R [ dxV^dx + R [ [f dxivg) + g dxiU)] dx = 
Jr ^ Jr ^ Jr 

for all (e, r?) G 



(2.30) 
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For X G M, we have 
(x) - E{dx(p{x)) - dxE{dx(f{x)) {x - dx^{x))\ -- 



<\\dll 



[x - y) d^E{y) dy 
{x - d^ip{x)f 



Multiplying the above inequality by foix), integrating over M, and using the definition of dx^p yield 

WiifJo) 



■(x) - E{dxipix)) - dxE{dx(p{x)) {x - dxi^ix))] fo{x) dx 



2.31) 



Owing to (2.30) with (^,?/) = {dxE,Q) and the property / = dxf#fo, we deduce that 



- [ if-fo) Edx-{1 + R) [ l^dlEdx-R I gdxUdxE) 



dx 



-2"^'""°° r 



Taking into account that {f,g) G i7^(M;IR^) by Lemma 2.1, we arrive, after integrating by parts once, 
to (2.15). A similar argument leads to (2.16). □ 

We next develop further an argument from the proof of [15, Proposition 2] which allows us to gain 
regularity on / and g by using the Euler-Lagrange equation. 

Corollary 2.3. The functions a// dx{{l + R)f + Rg) and ^ dx{f + g) both belong to L^(M) and 

dx[{l + R)f + Rg] < W2{fJo), (2.32a) 



r R^\\^dxif + 



< W2ig,go)- 



(2.32b) 



It is worth mentioning here that the estimates (2.32) match exactly the regularity of {f,g) given by 
the dissipation in the energy inequality (1.3). 

Proof. Consider G C^(M). We infer from (2.30) with = that, after integrating by parts, 
/ [(1 + R) f dxf + Rf dxg] idx = - [{x- dxifix)) (e o dxy^){x) /o(x) dx . 

jR T JR 

Since / = dx(pi^fo, it follows from the Cauchy-Schwarz inequality and (2.17a) that 

{x - dx<f{x)) (C o dx(p){x) fo{x) dx 



< 



(x - dxipix))"^ /o(x) dx 



1/2 



{^odx(pf{x) fo{x)dx 



1/2 



<W2{fJo) ( [ 

\Jr 



^^(x) /(x) dx 



1/2 



Therefore, 



/ [{I + R) f dxf + R f dxg] idx 

JR 



< 



W2UJ0) 



fix) f{x)dx 



1/2 



(2.33) 



Consider next a nonnegative function x ^ C'q°(]R) such that ||x||i = 1 ^^nd define Xni{x) := mximx) for 
m > 1 and x G M. Then, {Xm)m>i is a sequence of mollifiers in R and, given G Cq°{M) and m > 1, 
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the function t?/(m-^/^ + Xm* fY''^ belongs to C^(M). Taking ^ = t?/(m-^/^ + Xm* fY^"^ in (2.33), we 
obtain 



f d^[{l + R)f + Rg\ 
\lmr^l^ + Xm* f 



•d dx 



< 



f 



1/2 



The previous inequahty being vahd for all G C^(M), a duality argument yields 



^ d^{(\-,R)f + Rg\ 



^ni-y^ + X,n*f 



< 



W2{fJo 



f 



^-1/4 f 



1/2 



(2.34) 



Now, since / E H^(R) by Lemma 2.1, we have ||xm * / — /||oo < C'x lli^aj/lb m ^^"^ for some constant 
> depending only on x from which we deduce that 



/ 



m-i/4 + * / 
In particular, for x E M, 



< 



f -Xm* f 



^-1/4 f 
fix) 



Xm* f 



< 1 + 



and 



so that 



fix) 



j^-l/4 f 

< (i + II9../1I2) \/7M E L 

{ = v7ix) if /(x) = 0, 



m 



1/4 



(2.35) 



lim 

m— >oo 



Vm-^^^ + Xm * fix) 

f 



[ if fix)>0, 



f in 



\/m-V4 + x™ * / 

by the Lebesgue dominated convergence theorem. Since {1 + R)f + Rg belongs to H^{R) by Lemma 2.1, 



we conclude that 



/ 



d,[{l + R)f + Rg]^y^d,[{l + R)f + Rg] in L\ 



(2.36) 



Owing to (2.35) and (2.36), we may let m — > 00 in (2.34) and deduce that ^/f dx[{l + R) f + Rg] E 
and satisfies (2.32a). The proof of (2.32b) is similar. □ 

2.3. Interpolation. Thanks to the results established in the previous sections, we are now in a position 
to set up a variational scheme to approximate the solution to (1.1). More precisely, given (/o,5'o) £ ^2 
and T E (0, 1), we define inductively a sequence (/",5")„>o as follows: 



(/0,<7?) := (/o,<7o), 
(/;+!, 5^^) := inf T:^{u,v), 



with 



^iu,v) :-- 



2r 



(u,j;)e/C2 



Wi iu, /;) + Wi {v, g"^) ) + £iu, v) , (u, v) E IC2 



(2.37) 
(2.38) 
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the existence and uniqueness of (/""''^, (7""^^) being guaranteed by Lemma 2.1 for each n > 0. We next 
define two interpolation functions fr and Qr by (1.10), i.e. /r(i) := fr ^^'^ 5r (*) := 5™ * ^ ['^'^^ (n+l)r) 
and n > 0. By Lemma 2.2, we have 



JR JR 

[ [g^-gr^) Cdx + rR^ [ g^ d,{g!^ + g:^) d^^dx 
Jr Jr 



< 



zSlloo ttt2/ fn fn-1 



(2.39) 



for all n > 1 and G C^(M). Given T > arbitrary, we set := [T/t]. Summing both equations of 
(2.39) from n = 1 to n = A^, we find 

TV 



r i-{N+1)t f- 

/ {fr{T) -fo)Cdx+ fr d, {{1+R)fr + R <7r) dxdt 

JR Jt jr 

I- /.(7V+l)r p 

/ {9t{T) - go) ^dx + / gr [fr + gr) dxidxdt 

JR Jt jr 



< 



\m\\c 



n=l 



N 



(2.40) 



Y.wiigr,gr'] 



n=l 



(2.41) 



3. Convergence 



We gather in the next lemma various properties of the interpolations {frjgr) defined in Section 2.3 
which are consequences of Lemma 2.1 and Corollary 2.3. 

Lemma 3.1. There exists a positive constant Ci depending only on R, R^, /q, and go such that, for all 
t >0 and T G (0, 1), we have 



(n) 

(Hi) 
(iv) 

(v) 
(vi) 
{vii) 



fr{t) dx = gr{t) dx = 1, 

J] [w^K/;, /;-') + w^K^;^, ^r^)] <Cir, 

n=l 

£{fr{t),gr{t)) < £{fr{s),gr{s)), S G [0,i], 
/ {fr + gr){t,x) X^dx<Cl (1 + t), 

Jr 

^ [\\d.fr{s)\\l + \\d.,gr{s)\\f\ ds<Ci (1 + t), 

poo f 

/ / fr \d, [(1 + R) fr + R grf dxds < Cl, 
Jt Jr 

/ gr \dx {fr + gr)] dxds < Cl . 

Jt Jr 



(3.1) 

(3.2) 

(3.3) 
(3.4) 

(3.5) 
(3.6) 
(3.7) 



Proof. The property (3.1) readily follows from the fact that {fj^,g^) G IC2 for all n > and r > 0. Next, 
for r > and n > 1, the minimizing property of (/", (7") ensures that 



£{fr.9r)+ 



2t 



Wi{f:,fr') + ^Wiig!^,gr'' 
n,j 



<£{fr\gr'' 



{31 
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Given t G (0, oo) and s E [0,t\, we set N := [t/r], v := [s/r], and sum (3.8) from n = z^ + lupton = A^ 
to obtain, since {f^,g^){t) = {f? ,9^) and {fr,gT){s) = (/r,ffr)> 

N 



S{fAt),gAt)) + ^ ^ 



n=u+l 



<£{fr{s),gAs)). 



(3.9) 



The monotonicity property (3.3) is a straightforward consequence of (3.9) while the nonnegativity of £ 
and (3.9) with s = u = give 

N 



E 

n=l 



R 



<2£{fo,go) T. 



Since the right-hand side of the above inequahty does not depend on iV, we obtain (3.2). In order to 
prove (3.4), we combine (2.8) and (3.2) and obtain for t > with N := {t/r] 

frit, x) x^dx = I f?{x) x^dx<2 [ fo{x) x^ dx + 2T^|(/f , /o) 

N 

< 2 / Mx) x^dx + 2N ^W^if^jr') 

•J^ n=l 

[ fo{x) x^ dx + 4Nt £{fo, go) <C{l + t). 



< 2 



We next infer from (2.3) that, for n > 1, 



R 



r (||a,/;||2 + i?||a.(/; + 5;^)||2)<^(/-i)_^(/;) + _ {H{grr')-H{g^^)). 
Let N > 1. Summation from n = 1 to iV yields 

.{N+1)t jy 

(l|5x/r(s)||i + R \\d:c{fr+gr){s)\\l) ds < H{fo) - H{U{Nr)) + — {H{go) - HigriNr))) . 

(3.10) 

It now follows from Lemma A.l, (3.1), (3.4), and (3.10) that 

[■{N+l)r 

J {\\d,fris)f2 + R\mfr+gr)is)\\l) ds 



< ^(/o) + ^ Higo) + i^±^ 



+ ^(1 + x") lyfriNr) + griNr)^ < C (1 + iVr) , 



which entails the vahdity of (3.5) for t G [iVr, (A^ + l)r). 

We finish the proof by showing (3.6) and (3.7). By Corollary 2.3, we have for n>l 

//^ [(1 + R) /; + R < w^Ur. /r') • 

Summing over n > 1 and using (3.2) give 



Y,r^ ^/J^dA{i + R) f^ + Rg"^] <J2^2ifrJr')<CiT, 



n=l 



n=l 



whence (3.6). The proof of (3.7) also relies on Corollary 2.3 and is similar. 



□ 
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3.1. Compactness. We now turn to the compactness properties of (/^)^>o and [g^)^^Q and point out 
that the nonhnearity of (1.1a) requires strong compactness. We first observe that the compactness with 
respect to the space variable x is granted by (3.5) thanks to the following lemma. 

Lemma 3.2. The spaces F^(M)nLi(M, (l+x2)(ix) andL'^{M.)r\L^{M.,{l+x'^)dx) are compactly embedded 
in L^(M) and ff~^(M), respectively. 

Proof. Let us first consider a bounded sequence (/ij)j>i in //^(M) Pi L-'^(IR, (1 + x"^) dx). On the one hand, 
since H^iW) is continuously embedded in L°°(M) and C'^^'^iM), the Arzela-Ascoli theorem implies that 
there are h G i?^(M) and a subsequence of (/ij)j>i (not relabeled), such that (/ij)i>i converges to h in 
C([— i?, i?]) for all i? > 0. On the other hand, using once more the embedding of H^(W) in L°°(]R), we 
have for R> 1 

I \hi{x) — h{x)\^ dx < I \hi{x) — h{x)\^ dx + j \hi{x) — h{x)\'^ dx 

Jr J{\x\<r} j{|x|>_r} 

<2R \\hi - + ^ ll^i - ^lloo / x"^ \hi{x) - h{x)\ dx 



<2R - /i|lc([-i?,/?]) + sup|||/ii||oo 



x^ \hi{x) \ dx 



Letting first i — t- oo and then i? — t- oo shows that (/ii).(>i converges to h in L^( 

Next, let (/jj)j>i be a bounded sequence in L2(M)nLi(]R, (l + x^) dx) and denote the Four ier transform 
of hi by Thi for i > 1. A straightforward consequence of the bounds for (/ij)j>i is that (-F/ii)i>i is 
bounded in L^(M) H W'^'°°{W). Arguing as above, this implies that (J^/ij)j>i is relatively compact in 
L^(]R, (1 + x'^)~^ dx). Coming back to the original variable, (/ii)i>i is relatively compact in if~^(IR) as 
claimed. □ 



We next turn to the compactness in time and prove the following result: 

Lemma 3.3. There is a positive constant C2 depending only on R, R^, /q, and qq such that, forr G (0, 1) 
and {t,s) £ [0, 00) x [0,oo), 



Wfrit) - /r(s)|U-3 + Ibr(i) " <7r(s)|k-3 < C2 V\t-s\+T. 



(3.11) 



Proof. Consider t G (0,oo), s G [0,t], and define the integers N := [t/r] and v := [s/t]. Either N = v 
and fr{t) - fr{s) = satisfies (3.11) or iV > i/ + 1 and it follows from (2.39) that, forn G {z^ + 1, • • • , N] 
and i G C^(M), 



r /■(n+l)r i- 

/ (/; - /;-') idx < / fr{s) \d. [(1 + R) f, + R g,] {s)\ \d^i\ dxds 

Jr Jut Jr 



+ \Mlh^wi{f?,fr') 
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Summing the above inequality from n = i' + l to n = N and using (3.1), (3.3), (3.6), and the Cauchy- 
Schwartz inequahty, we are led to 



(frit) -Ms)) Cdx 



N 

^ E 



[ (/; - fr') ^ 

JR 



dx 



+ 



AN+1)t r 

< / / fr{s) \d, [(1 + R) fr+R gr] {s)\ \d^^\ dx ds 

J{u+1)t Jr 



< 



\\dM\c 



N 

n=v+l 
{N+l)r 

\\fr{s)\\T 



Tr [(1 + R) fr + R gr] (s) 



ds 



{Vt 



T + T 



< C U\\w2,o. [Vt-S + T + T) 

Since H'^(W) is continuously embedded in VK^'°°(R), the claimed estimate for fr{t) — fris) follows by a 
density argument. A similar computation relying on (2.39), (3.1), (3.3), and (3.7) gives the same estimate 
for gr{t) — gr{s) and completes the proof of Lemma 3.3. □ 

We are now in a position to establish the strong compactness of [fr, gr)r>o hi L^((0,r) x M) for all 
T > as announced in (1.11). 

Lemma 3.4. There are a sequence {Ti^)k>i, Tk ~^ 0, and functions f and g in C {[0, oo); (JR.)) such 
that, for all t > 0, 

{f,^{t),gr,{t))^{f{t),g{t)) m H-'\R;M.^) , (3.12) 

(/r,,5rj ^ (/,5) L\{0,t)xR;R^), (3.13) 

{fr,,gr,)^{f,g) a.e. m (0,oo)x]R. (3.14) 

Proof. The proof relies on [3, Proposition 3.3.1] and [17, Lemma 9]. Indeed, it first follows from (3.1), 
(3.3), (3.4), and Lemma 3.2 that {fr{'t))r&{o,i) li^s in a compact subset of H~^(M.). This fact, together 
with Lemma 3.3 and a refined version of the Arzela-Ascoli theorem [3, Proposition 3.3.1] ensures that 
there are a sequence (rfc)fc>i, — )■ 0, and a function / G C([0, oo); i?~^(M)) such that (fr^it)) converges 
towards f(t) in i7~'^(M;M^) for each t > 0. Since the same argument applies for (5T)rG{nii' we have 



in H- 



=(0,1)' 

the convergence (3.12), 



established (3.12). We then infer from (3.3), the embedding of 
and the Lebesgue dominated convergence theorem that 

{fr„gr,)^{f,g) in L\0,T;H-^{R;R^)) for ah r>0. (3.15) 

Now, given 6 G (0, 1) and T > 1, the estimates (3.1), (3.3) (with s = 0), (3.4), and (3.5) in Lemma 3.1 
ensure that 

{frk,gTk)k>i is bounded in L'^{6,T; H^{R) n L^{R, [1 + x"^) dx)) . (3.16) 

Since i?^(M)nL^(]R, (1+x^) dx) is compactly embedded in L^(M) by Lemma 3.2 and L^(M) is continuously 
embedded in iJ~^(R), we are in a position to apply [17, Lemma 9] and deduce from (3.15) and (3.16) that 
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{frk,9Tk)k>i converges towards {f,g) in L^{{6,T) 



^). Owing to (3.3), this convergence may actually 



be improved to (3.13). The a.e. convergence (3.14) then follows from (3.13) after possibly extracting a 
further subsequence. □ 

Finally, (3.5) implies that, after possibly extracting a further subsequence, we may assume that 

{d,fr„d,gr,)^{d,f,d,g) in L\{6,T)xR) for all 0<6<T. (3.17) 
Now, combining (3.6), (3.7), (3.13), and (3.17), we obtain 

a,, [(1 + ii) + ii <7rJ - y^d^il + R) f + Rg] 



Vff^ 9x (/rfc + 5rJ ^ Vg d^{f + g) 

for < (5 < T, while (3.13) and (3.17) imply that, for < 5 < T, 



in L\{5,T)x 



(3.18) 



[fr.dA{l + R)fr.+R9r.] ^ f d [(I + R) f + R g] .^^ ^1 ^T) X M) . (3.19) 



3.2. Passing to the limit. Combining the convergence (3.12) with the estimates (3.1), (3.3) (with 
s = 0) and (3.4) in Lemma 3.1 ensures that {f{t),g{t)) G IC2 for all t > 0. Moreover, gathering (3.3), 
(3.5), (3.13), and (3.17), we conclude that {f,g) satisfies the integrability properties (i) of Theorem 1.1. 
In addition, it follows from (3.12) and Lemma 3.3 that 



Wfit)- f{s)\\H-^ + \\9{t)-9{s)\\H-s <C2 V\t-s\, (t,s)G [0,00) X [0,00), 



(3.20) 



which proves assertion (ii) of Theorem 1.1. 

In order to establish the estimate (b) of Theorem 1.1, we pick T > and set A''^ := [T/r^] for all 
integers k > 1. Then, we infer from Corollary 2.3 and (3.9) (with s = 0) that for all /c > 1 we have 



1 



'Tfc 

n=l 



fr,{a)dAil + R)fr,+R9r,]i<T 



M^|(/-,/--i) R Wi{gl,g!^-^' 



2rfc 



+ 



Rn 



2Tk 



+ RRf, y 9Tk (c^) [frk + 9Tk ] (f^) 
<£{fo,go)-£ifr,iT),gr,{T)). 



da 



Letting k — )• 00, the convergences (3.13) and (3.18) lead us to 



+ RRii 



Vgia) d,[f + g]ia) 



y/J{a) [(1 + R)f + Rg] (a) 
<£{h,go)-£{f{T),g{T)). 

for all 5 G (0, 1), whence the desired assertion (5) of Theorem 1.1 after letting — )■ 0. 



da 
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Now, we identify the equations solved by / and g. To this end, fix E C^(M), t E (0,oo), s E (0,t) 
and set N := [t/r] and u := [s/t]. We infer from (2.40), (3.1), (3.2), and (3.6) that 

/ (frit) - fris)) idx+ f f U{a) [(1 +R) fr + R Qr] {^) dxda 

Jm. Js Jr 

r r{N+l)T r 

/ (frit) -fo)^dx+ / fria) [(1 + R) U + R gr\ (t) d^^ dxda 

Jm. Jt Jr 

/ (fris) - fo) idx + /^'^'^" / fr{a) [(1 + R) U + R gr] {(y) d.^idxda 
Jr Jt Jr 

f{N+l)T f- 

/ / fr{a) [(1 + R) U + R gr] {t) d^i dxda 

Jt Jr 

[ U{(y) dx [(1 + R) U + R 9r\ {<y) d^C dxda 
Js Jr 



+ 



+ 



N 



<\\diaoo ^i^i(/;,/r') 



n=l 



+ \\dM\oo [^'^^'^^ \\U{^)\\Y^ yq;dA{i + R) U + Rgr]{c7) 

J s 



+ m\\c 



/■^^+^)" 1/2 ^ 

J ll/r(^)llf ^/U^A{l + R) fr + R9r]{<T: 



da 



da 



Taking r = in the above inequahty and letting /c — )• oo with the help of (3.12) and (3.19), we end up 
with the first identity in (1.12) 

/ (fit) - f{s)) ^dx+ f [ f{a) [{1 + R) f + R g] (a) d^i dxda = . 
Jr Js Jr 

The proof of the second one being similar, it remains to check the property (a) stated in Theorem 1.1. 
To this end, we first claim that 



(^,ln^,,<7,,ln<7,J ^ (/ln/,5ln5) in L1((0, T) x M) , T>0. 
Indeed, by (3.13) and the continuity of r i— t- r In r in [0, oo), we have for T > 

(/rfc In frk , ark In fi-r^ ) — > {flnf,g In g) a.e. in (0, T) x M. 
Moreover, it readily follows from (3.3) (with s = 0) that 



(3.21) 



(3.22) 



{Uk In Uk > 9Tk In ^Tfe )fc>i is uniformly integrable in ((0, T) x R; M^) , 



(3.23) 
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while (3.3), (3.4), and the inequahty |rlnr| < max{r, 1}, r > 0, guarantee that, for R> 1, 
(■T f f-T 



[ [ \fr,lnUJdxdt < 2 [ [ yq7, 1[0,1]{ fr,)dxdt 

Jo J{\x\>R} Jo J{\x\>R} 

Jo JUx\>R\ 



10 J{\x\>R} 
t-T 



+ '^11 fU^l[R,oo){fr,)dxdt 
Jo J{\x\>R} 



T 



1/2 / „^ „ \ 1/2 

dxdt \ 



< 2 / / X Vrfe dxdt I / / 

\J0 J{\x\>R} J \J0 J{\x\>R] X 

+ 2VR [ [ fr, dxdt + ^ / / /2 dxdt 

Jo J{\x\>R} VR Jo J{\x\>R} 

<cl±l. (3.24, 



Due to (3.22)-(3.24), we are in a position to apply Vitali's convergence theorem (see, e.g., [11, Theo- 
rem 2.24] or [13, Theoreme 1.4.13]) and deduce the claim (3.21) for {fTk)k>i: the proof for ((7rfc)fc>i being 
identical. Consequently, after possibly extracting a subsequence, we have also 

{H{U,{t)),Higr,m ^ iHifit)),H{gm a.e. in (0, oo) , (3.25) 

the functional H being defined in (2.4). We next infer from (3.17) and the Fatou lemma that, for t > 0, 

[\\\dxf{s)g + R\\dM + 9)is)\\l) ds=hm [\\\d,f{s)\\l + R\\dM + 9)is)\\l) ds 
Jo o->0 Js 

<lim hminf [\\\d,fr,{s)\\l + R \\dMr, + 9r,){s)g) ds . (3.26) 

Owing to (3.25) and (3.26), we may pass to the limit as A; — )• oo in (3.10) to obtain the assertion (a) of 
Theorem 1.1, which completes its proof. 
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Appendix A. Some technical results 
We first collect some well-known properties of the functional H defined by (2.4). 
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Lemma A.l. Let h be a nonnegative function in (1 + x'^)dx) n Then hlnh G L^(M) and 

there is a positive constant Cg such that 

(A.l) 
(A.2) 



h{x) \ In h{x)\dx < Ci+ / h{x) (l + x^) dx + \\h\\l , 

H{h) > -Ce- [ h{x) (l + x^) dx . 
Jr 



Proof. Introducing the function u}{x) := e (^+^^), x G M, and using the monotonicity of r i— ?• r| lnr| in 
[0, 1/e], we have 



h{x) \ In h{x) \ dx 



{h{x)<uj{x)} 



h{x) \ In h{x) \ dx + 



{uj{x)<h(x)<l} 



h{x) \ In h{x) \ dx 



+ 



{h{x)>l} 



h{x) \ lnh{x) \ dx 



< 



{h{x)<uj{x)} 



(l + x- 



')dx+ / 

J {uj{x 



)<h{x)<l} 



h{x){l +x'^)dx 



+ 



{h{x)>l} 



h{x){h{x) — 1) dx 



< I e-(i+NP) (i^x^)dx+ [ h{x)il + x"^) dx + 
Jr Jr 



whence (A.l). Similarly, 

H{h) > [ 
Jh 



{h{x)<uj{x)} 



h{x)\nh{x) dx + I h{x)\nh{x) dx 

J {u}{x)<h{x)<l} 



> 



I {h{x)<u)(x)} 

from which (A.2) readily follows. 



I e-(^+l^l') (l + x2)(ix- /" h{x){l + x"^) dx , 

J \h{x)<uj{x)\ J {'^(x)<h(x)<l} 



□ 



The next results allowed us to identify the limit of some terms arising in the derivation of the Euler- 
Lagrange equation in Lemma 2.2. 

Lemma A.2. Consider h G i7^(M) and ( G C^(M). Setting Q := id+e ( for e > 0, we have 

hoQ — >h in L^m and hoCe-h _^ ^^^^ L^m. (A.3) 

Proof. Since C,^ is a C°°— diffeomorphism from M onto M for e small enough, its inverse C,~^ is well-defined 
and satisfies 

\x-Q\x)\<e\\CU, x^W. (A.4) 
It follows from the Cauchy-Schwarz inequality, the Fubini theorem, and (A.4) that 

2 



\\hoCe-h\\l = [ I ['^^^'^^ dMy)dy] dx < [ \x - Ce{x)\ [^^^^\dMy)\^ dy 

Jr \Jx J Jr Jx 



< e 



\dMy)? \y-Q\y)\ dy<e' IICIlL Wml, 



dx 
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dxh{y)dy C{x) — ldxh{x) C(x) 



which gives the first assertion in (A. 3) and the boundedness of {{ho Q^, — h)/e)^ in L^(M). Next, since 
h G i/-^(]R), almost every x G M is a Lebesgue point for dxh and, for such points, 

h{x + eC{x)) - h{x) _ 1 r^+ea^) 

Therefore, {{h o Q — h)/e)^ converges a.e. to C,dxh as e — )• and is bounded in and the second 

assertion in (A. 3) readily follows from these two facts. □ 

The first assertion of Lemma A. 2 is actually true in a more general setting: 

Lemma A. 3. Consider h G H^{W) and a sequence (Ce)e>o of functions in C^(M) such that ujs '■= 
\\Ce ~ id lloo — > as e — )• 0. Then 



hoQ — >h m L2(M). 

e-i-0 

Proof. As in the proof of Lemma A. 2, it follows from the Cauchy-Schwarz inequality and Fubini's theorem 
that 

\dxh{y)\'^dy dx < LVe / / \dxh{y)\'^ dydx 



hoCe-h\\l < I \x-Ux)\ 



< 2 \\dxh\\l , 

and the right-hand side of the above inequality converges to zero as e — )■ 0. □ 
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